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ABSTRACT 

Our heuristic understanding of the abundance of dark matter halos centers around the concept of a density 
threshold, or "barrier", for gravitational collapse. If one adopts the ansatz that regions of the linearly evolved 
density field smoothed on mass scale M with an overdensity that exceeds the barrier will undergo gravitational 
collapse into halos of mass M, the corresponding abundance of such halos can be estimated simply as a fraction 
of the mass density satisfying the collapse criterion divided by the mass M. The key ingredient of this ansatz 
is therefore the functional form of the collapse barrier as a function of mass M or, equivalently, of the variance 
a^{M). Several such barriers based on the spherical, Zel'dovich, and ellipsoidal collapse models have been 
extensively discussed. Using large-scale cosmological simulations, we show that the relation between the 
linear overdensity and the mass variance for regions that collapse to form halos by the present epoch resembles 
expectations from dynamical models of ellipsoidal collapse. However, we also show that using such a collapse 
barrier with the excursion set ansatz predicts a halo mass function inconsistent with that measured directly in 
cosmological simulations. This inconsistency demonstrates a failure of the excursion set ansatz as a physical 
model for halo collapse. We discuss implications of our results for understanding the collapse epoch for halos 
as a function of mass, and avenues for improving consistency between analytical models for the collapse epoch 
and the results of cosmological simulations. 

Subject headings: 



1. INTRODUCTION 

A central concept in the modern theory of galaxy for- 
mation is the connection between characteristics of the lin- 
ear density field and the abundance and properties of viri- 
alized dark matter halos in the contemporary universe. The 
power spectrum of density perturbations seeded by infla- 
tion (e.g., Guth & Pi 1982; Bardeenet al. 1983; Stai-obinskv 
1983) and the cosmological transfer function (e. g., Peebles 
1982: .Bardeen et al.ll 19861 lEisenstein & Hull 19981) determine 



the character of the subsequent nonU near growth of struc - 
ture through gravitational clustering (jWhite & Rees' Il978l) . 
Growing perturbations in the i nitial den sity field ser ve as the 
sites of galaxy formation (e .g., lPeeblesllT965l: ISachs & Wolfd 
119671 IWhite & Reel 119781) . In the context of a cold dark 
matter cosmology, these processes give r ise to the char- 
acteri stic mass scale of observed galaxies dRees & Ostrikerl 
119771: 'Blumenthal et al. 1984). Cosmological observations 
have both motivated and verified this picture, most recently 
with measurements of galaxy clustering (e.g., i Percival et al.l 
120071) . the linear power spectru m of cosmological struc- 
tures (e.g, iMcDonald et al.ll2006h . and high-precision mea- 
surements of the cosmological microwave background radi- 
ation (e.g., Dunkley et al. 2008). 

Methods for calculating the abundances of nonlinear, col- 
lapsed structures have been developed to link the growth 
of density perturbations with the observed number densi- 
ties of galaxy- and cluster-scale objects. Dynamical models 
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for the collapse of individual dense patches into viri alized 
structures, such as the spherical collapse (Gunn & Go t ljl 19721 
see Appendix lAli and ellipsoidal collapse (Eisenst ein & Loebl 
1995; Bond & Myers 1996, see AppendixlBt models, provide 
physically-motivated methods for estimating the necessary, 
linearly-extrapolated overdensity (the "collapse barrier") for 
a region to break from the cosmic expansion, condense, and 
form a high-density, virialized structure (i.e., a dark matter 
halo). When combined with the statistics of the initial den- 
sity field, the collapse barrier can thereby be utilized to es- 
timate the abundance of dark matter halos as a function of 
mass and redshift. The purpose of this paper is to re-examine 
the connection between the collapse barrier and halo abun- 
dance, and test the common assumptions and methodologies 
used to calculate the mass function of dark matter halos from 
a dynamic al model for their c ollapse (e.g., the "excursion set" 
f ormahsm,|Bond et al.llT99l1) . 

iPress & SchechteJ ('1974') first used the spherical collapse 
model to calculate the abundance of galaxies. They assumed 
the probability distribution function dP((5R)/d(5R (PDF) of the 
smoothed overdensity field (5r = {pm- Pm) / Pm, where the mean 
matter density is p,„ and the density is averaged over a region 
of typical size R containing mass M cx pmR^, was a Gaus- 
sian with a scale-dependent variance cr^iM). They integrated 
this Gaussian PDF above the typical collapse overdensity 
5c (i.e., P{5m > Sc) (x erfc[-(5c/V2(T(M)]) and differentiated 
with respect to mass M (i.e., dP/dM cx exp[-5^/2cr^(M)] x 
dCT"'/dM) to arrive at th e fraction of all mass co ntained in 
objects of mass M. The Press & SchechteJ 11974) calcula- 
tion accounts for only half of the total universal mass den- 
sity in bound objects because it does not address underdense 
regions contained within still larger regions for which the 
threshold 5m > Sc is satisfied. To remedy this shortcoming, 
iPress & SchechteJ ([1974) multiplied their final answer by a 
fa ctor of two t o acco unt for all mass with little justification. 

iBond et aH (Il991l) studied the properties of sets of regions 
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above the threshold, the "excursio n sets" of the dens ity field. 
Using the excursion set formalism. lBond et alJ (1199 ll) demon- 
strated that by filtering the initial overdensity field on avari- 
ety of mass scales, the results of iPress & SchechleJ d 19741) 
could be derived in a manner that accounts for patches of low 
density embedded in large, high-density regions collapsing 
on larger s c ales ( the " cloud-in-cloud" prob lem). Moreover, 
iBond et alj (1199 1') and'Lacev & Cold d 19931) showed that the 
excursion set formalism provided a means to compute other 
halo properties such as their mass acquisition histories. The 
ex cursion set theory of halo abundance was later extended 
by iMo & White! (1 19961) to describe their spatial clustering 
through the bias parameter = f /,/, /^,„ relating the halo {£^i,h) 
and mass correlation functions (a computation that re- 
covers the "peak-background split" result developed in lKaised 
[1984; Efstathiou et al, 1988; Cole & Kaiser 1989^). The de- 
tails of t he excursion se t theory are collected in the recent 
review bv' Zentnerf (120071) . 

Numerical simulations of cosmological structure forma- 
tion, which were developed concurrently with the analyti- 
cal collapse calculations, demonstrated that nonlinear gravi- 
tational collapse produces hal o mass functions that are incon- 
sistent with the predictions of Press & Schechted(ll974t) . Evi- 
dence for this disagreement, as well as corresponding discord 
in the spatial clustering of halos, developed over twent y years 
(ej.,lEfstathiou et al. 1988; White et al. 1993; Lacev & Cold 
1994; Eke et al. 1996; Gross et al. 1998; Tormen 1998; Jing 
1998; Lee & Shandar ia 1999; Jing..! 999; .Porciani et al..,1999; 
Governato et al. 19991). The inability of the simple, spherical 
collapse model set within the excursion set formalism to de- 
scribe the results of cosmological simulations and the need for 
robust predictions of the abundance of halos for comparisons 
with observations led to accurate formulae for mass func- 
tions determined by numerical fits to the results of N-body 
simulations (e.g..'Sh eth & Tormenlll999Hjenkins et al."200l'; 
Warren et al. 2006, for more recent results see 0"inker et al. 
20081) . These fitting formulae were not the results of spe- 
cific dynamical models for the collapse of dark matter ha- 
los. Rather, they were developed through a practical approach 
of trying to reproduce accurately the results of cosmological 
simulations. 

Additional dynamical models for the growth of structure 
were developed in an attempt to better reproduce simu- 
lated mass functions while retaini ng a comparably s imple , 
physically-motivated framework. iLee & ShandarinI (1 19981) 
presented r esults for the halo mass function motivated by the 
IZeFdoviciil Il97 0) pancake collapse model, which they ex- 
tended to account fo r the col l apse o f halos along each of their 
principal axes. Sh ethetal.1 dioOlL SMTOl) ar gued that the 
functional form of the ISheth & TormenI (I1999L ST99) mass 
function can be motivated by the ellipsoidal collapse model of 
iBond & MversI (!l996). For a given ellipticity e and prolaticity 
p of th e shear field about an overdensity S, the lBond & MversI 
dl996l) model provides a method for estimating the linearly- 
extrapolated overdensity at collapse (the "ellip soidal collapse 
barrier "). SMTOl built on a calculation by iDoroshkevichl 
(11970) to determine a probability distribution of shear ellip- 
ticities and prolaticities, which was then used to find the most- 
probable values for e and p as a function of a(M). Combining 
these results, SMTOl found an effective ellipsoidal collapse 
barrier as a function of the "peak height", = Sc/a-{M), alone. 
With this new, mass-dependent collapse barrier SMTOl com- 
puted a halo mass functi on using the excursion set prescrip- 
tion of IBond et al.l dl991h and showed that the predicted func- 



tional form was close to that measured in cosmological sim- 
ulations by ST99. The ST99 mass function has therefore be- 
come associated with the ellipsoidal collapse model. 

Note, however, that the ellipsoidal collapse model predicts 
that the collapse barrier converges to the spherical collapse 
barrier Sc = 1.69 (for il„ = 1), in the high-mass limit, be- 
cause the rarest peaks have a preferentially spherical shape. 
SMTOl, on the other hand, found that the collapse barrier had 
to be lowered in the high-mass limit to ^/osrSc ^ 0.84(5c, in or- 
der to reproduce the mass function m easured in cosmological 
simulations (see also Shet h & T ormen 2002). Although they 
argued that the lower value can be motivated by the mass defi- 
nition of the Friends-of-Friends algorithm which they used to 
identify halos in simulations, this rescaling is not well justi- 
fied. Indeed, it disagrees with the fact that the collapse of the 
highest mass halos should be described well by the spherical 
collapse model. We argue below the lowering of the collapse 
barrier is instead required by the internal inconsistencies of 
the excursion set ansatz and explicitly demonstrate that such 
inconsistencies persist for a wide variety of halo mass defi- 
nitions, including the Friends-of-Friends and spherical over- 
density criteria (see Appendix iDl for a detailed discussion). 

In the context of the excursion set formalism, the adoption 
of a shape for the collapse barrier (along with choice of pre- 
scription with which to smooth the density field) effectively 
determines the mass function. In this paper we test the ex- 
cursion set ansatz by measuring the effective collapse barrier 
for halos formed in cosmological simulations and comparing 
excursion set predictions for the halo abundance given such 
a barrier to the halo mass functions measured in the simula- 
tions. 

The organization of this paper is as follows. In ^ we 
review the theoretical background of the excursion set for- 
malism, common dynamical models, and collapse barriers 
from the literature. In ^ we measure the linear overden- 
sity of collapsed regions in cosmological simulations to study 
the consistency between the simulated collapse of halos and 
the predictions of dynamical models. We then compare the 
abundance of dark matter halos predicted from the excursion 
set formalism and dynamical models with the simulated halo 
mass function in ^ We discuss our results in ^and summa- 
rize our results and conclusions in ^ The paper contains four 
appendices: a review of the s pherical col l apse m odel (Ap- 
pendix [A]), a review of the Bo nd & MyersI dl996|) ellipsoidal 
collap se model (AppendixlBl). a summary of the'Zhang & Hiiil 
d2006l) analytical method for calculating the excursion set 
mass function (Appendix ICll. and a study of the connection 
between the excursion set ansatz and the halo mass definition 
(Appendix iDb. Unless noted otherwise, throughout the paper 
we assume a flat A-Cold Dark Matter (ACDM) cosmology. 
Cosmological parameters in our calculations and simulations 
are close to the values suggested by observations: Vl,„ w 0.3, 
f^A « 0.7, and Hq w 70 km s"' Mpc"' . 

2. THEORETICAL BACKGROUND 

Consider a Gaussian random density field p(x) as a function 
of spatial location x with mean matter density p„,. At every lo- 
cation X we can define the overdensity (5(x) = [p(x) - p,„] / p,„ 
of the field and the overdensity smoothed on a comoving 
length scale as 

5fi„(x)= [dyS(x')W(\x-x'\,R^) (1) 
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where W(x,/?w) is a spherically-symmetric smoothing win- 
dow with a characteristic radius 7?w centered about loc ation x. 
Here, we follow the notation in the recent review by .Zentned 
(|2007) but use the symbol 6r„ to represent the overdensity 
field smoothed on a scale 7?w and reserve 6 to represent the un- 
smoothed overdensity field. The smoothed overdensity field 
Sr„(x) is also a Gaussian random field with a variance given 
by 

S(M) = aHM)={Sljx)) = j A\k)\W{k,R^)\^A\nk, (2) 

where A^(k) = k^P(k)/2TT^ is the dimensionless power spec- 
trum of density fluctuations with a wavenumber k and 
W{k, /?w) is the Fourier transform of the real-space smooth- 
ing window W(x,/?w)- The details of both the averaging pro- 
cedure and the mass-radius relation are fixed by the choice of 
filtering function used to smooth the density field. In CDM 
models the variance monotonically decreases with increasing 
length or mass scale. Consequently, once the filter function is 
specified this variance can also be used to label the size of the 
smoothing region, and we can write Sr^ as Sm or Ss where M 
is the mass contained within the window of length scale /?w 
and S is the variance a^(M). In what follows, we adopt the 
common practice of labeling the size of the smoothing region 
by either the length scale 7?w, the corresponding mass scale 
M, or the variance S, and use the labels interchangeably. 

If we consider a fixed location x and monitor the behavior 
of the smoothed overdensity (5s(x) as we decrease the mass 
smoothing scale M from some very large value (and hence 
increase the variance S from some value ^ S]), Ssix) will 
execute a (not necessarily Markovian) random walk where the 
smoot hed overdensity ^ ^(x) wil l satisfy the Langevin equation 
(e.g.. iBond etal.ll 19911 see also IChandrasekhaJ 1 943l) 

-^ = Q{\x,kW(k,R^), (3) 

where Q(\nk) is a Gaussian random variable with zero mean 
and variance 



AS 



(e^dn^)) = — - = A2(fe). 
^ ' dlnfc 



(4) 



The variation represents an ensemble of local realizations of 
the density field. 
For a sharp ^-space tophat filter of the form 



Wi(A;,«w) = 6(l-'t/?w), 

: jc<0 
e(jc) ={ \ : X = , 

1 : x>0 



(5) 
(6) 



the "trajectory" of overdensity as a function of variance (5s(x) 
recovered by integrating Equation[3]will be a Markovian ran- 
dom walk because each Fourier coefficient of the Fourier- 
transfo rmed density field is an independent random variable 
(e.g., Lacev & Cole 1993). More generally, the trajectory may 
vary more smoothly as the variance is increased if the Fourier 
transform of the real- space window function has broad side 
lobes in A:-space (e.g.. IZentnedl2007h . For instance, the real- 
space tophat filter, given by 



W,(R,R^) = ( ^R' ] 9(1 -R/R^), 



with Fourier transform 

WAk,Ry,): 



3 (sin kRw - kR^/ cos fc/?w) 
(kRw? ' 



(7) 



(8) 



has extended side lobes in fe-space that correlate the smoothed 
overdensities over a considerable range of smoothing scales. 
Hence, the integral of Equation [3] will vary more smoothly 
with the variance 5 if a real-space tophat is used rather than a 
^-space tophat. 

2.1. Collapse Barriers and Halo Formation 

The excursion set theory of halo abundance, clustering, and 
formation is based on an ansatz that the locations and sizes of 
virialized dark matter halos can be related to the properties of 
peaks in the initial density field at some very high initial red- 
shift Zinit ^ 1, when the density field is in the linear regime 
{5u <C 1). Almost universally, the specific form of this ansatz 
is that a region will collapse and form a dark matter halo if 
its smoothed overdensity, evolved forward in time from Zim. 
according to linear perturbation theory, exceeds some thresh- 
old value. Consider an object that collapses at some redshift 
Zc < Zinit- The linearly-extrapolated overdensity is 



SrJx,Zc) = SR^(x,Zink)D(Zc)/D(Ziaid, 

where D{z) is the linear growth function, given by 

(l+z')dz' 



D{z) = DqH{z) 



H\z') 



(9) 



(10) 



The Hubble parameter 



Hiz) = Ho [nui +zf+(i-n„-nA){i+zf+nAY^^ ( 1 1 ) 

describes the rate of change of the universal scale factor as 
H = a /a. The collapse condition is then simply 

SRjx,Znm)D(z,)/D(Zian) > B, (12) 

where B is referred to as the "collapse barrier" The excursion 
set ansatz is specifically that the largest smoothing scale /?w oc 
(M/p,„y/^ at any point x for which Equation fT2l is satisfied 
will collapse and form a halo of mass M at redshift Zc- 

The value of the collapse barrier B is usually determined by 
a dynamical model for the collapse of overdense patches in 
a background cosmological environment. In general, the col- 
lapse barrier need not be a single number and B can be a com- 
plicated function of the properties of the local linear density 
field (including its spatial derivatives), the smoothing window, 
and the smoothing scale. In the following sections, we study 
the spherical collapse and ellipsoidal collapse models for B. 

Having specified the collapse condition and the form of the 
collapse barrier, the number of collapsed objects at a given 
mass M or variance S will be determined by the probabil- 
ity distribution fiS)dS of variances where random realiza- 
tions of trajectories Ss(x), computed according to Equation 
[3] first cross the barrier B. This first-crossing distribution can 
be determined using a Monte Carlo procedure. By integrat- 
ing Equation [3] for many locations in the density field [cor- 
responding to many realizations of Q(lnk)], the first-crossing 
distribution may be approximated by a histogram of the bar- 
rier crossings as a function of S for the ensemble of trajecto- 
ries. The first-crossing distribution is often written as a func- 
tion of the peak height 



(13) 



where Sc = 1.686 is the linear overdensity for spherical col- 
lapse in an fl„, = 1 cosmology (see ^A] below for details). 
We will frequently change variables from S to i^c so that we 
can discuss and plot the first-crossing distribution in terms of 
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1 2 3 4 5 

ct(M)=i5 !v 

Fig. 1. — Some com mon collapse barrie rs discussed in the text: the 
spherical collapse barrier (Gunn & Gottf'WV?, 5c ~ 1.686 for f2,„ = 1, soHd 
gray line), the elHpsoidal collapse baiTier (Sheth et al. 2001, dashed red Hne), 
an alternative fit to the ellipsoidal collapse barrier ( Desjacqu es 2008, soUd 
purple hne), and the modified elhpsoid al collapse barrier associated with the 
ISheth & Toriiiei] fT999) mass function jSheth et ai]l200U dotted green hne). 
Note that while the ellipsoidal collapse barriers, as formulated, converge to 
sp herical collapse ba rrier at large masses [small (t(M)], the modified barrier 
of lSheth et alj i200TI) converges to a lower value of y/a^Sc ~ 0.845j.. 



from the collapse model is to assume that this equivalent lin- 
ear overdensity can be used to identify collapsed regions in 
the initial density field without needing to consider the field's 
nonlinear evolution. The model predicts that a region with 
initial physical overdensity ^^(zinit) will collapse when the 
linearly-extrapolated overdensity exceeds 

Ssiz) = Ss{z,ni,)D(z)/D(zMt) > Bsc = <5c. « 1 .686. (15) 

We refer t o the barrier 6e as the "spherical collapse barrier" 
lEke et al] (11996') provide analytical solutions for the value of 
Sc in cosmologies with 0„, + Ha = 1 . The spherical collapse 
b arrier is independe nt of mass scale and initial overdensity. 

iBond et aTl (Il991h used the excursion set formalism to cal- 
culate the first-crossing distribution associated with a barrier 
that is constant as a function of mass scale, such as the spher- 
ical collapse barrier. For a sharp-fc window function [Equa- 
tion ©I, the spherical collapse barrier first-crossing distribu- 
tion is 



2tt 



1/2 



exp 



(16) 



The function /sc is normalized such that J fsc(i^c)<ii'c = 1, im- 
plying that all mass in the universe is incorporated into col- 
lapsed objects. In this context, the normalization arises be- 
cause the variance S = cP-{M) increases monotonically toward 
infinity as M tends toward zero, while the barrier height 5c, 
remains fixed. As a result, any random walk will cross the 
barrier dr at some scale. 



The comoving abundance of halos (d«/dM) x AM in the 
mass range AM about mass M (the "mass function") is related 
to the first-crossing distribution f{vc) by 



AM = — fivc) 



dM 



AM. 



(14) 



This correspondence between the mass function and first- 
crossing distribution can be understood as a variable change 
from a distribution in peak height Vc to a distribution in mass 
M, with a normalization that accounts for the partitioning of 
mass elements into halos of mass M (i.e., p,„/M). The shape 
of the first-crossing distribution f{vc) is expected to be rela- 
tively independent of the cosmological model at low redshift 
because it depends primarily on the primordial power spec- 
trum, the collapse barrier B, and the smoothing window (e.g., 
ST99). The halo abundance dn/dM additionally depends on 
how the relation between the trajectory mass smoothing scale 
M and the actual halo mass is defined and the shape of the 
power spectrum P(k) (through dcr^/dM). When necessary, 
we will associate halos with roughly virialized regions of size 
^200 with a mean physical overdensity A = 200 (to ease com- 
parison with the mass definition used bv [Tinker et al.ll2008l 
which is defined relative to the background density, and not 
the critical density) and virial mass M200 = 47rAp,„/?^QQ/3. 
However, we explore variations on this definition for com- 
pleteness (see Appendix iDli. 

2.2. Spherical Collapse Barrier 

iGunn & Gotj (Il972h modeled the dynamical evolution of 
an overdense spherical region in a background cosmology 
(their calculation is detailed in Appendix |A|i. The collapse 
barrier is computed by evolving the density in the spherical 
region according to the linear theory until the time of collapse. 
The final step in calculating the excursion set halo abundance 



2.3. Ellipsoidal Collapse Barrier 

The spherical collapse model is likely too simplistic be- 
cause peaks in the linear density field ar e, in general, locally 
triaxial (see, e.g., iDoroshkevich 1970; .Bardeen et al.|[i986) . 
A number of dynamical collapse models designed to account 
for deviations from spherical symmetry have been explored 
(e.g., IZel'dovich 1970; Nariai & Fuiimoto 1972; Hoffma3 
'1986'; Bertschinger & Jain 1994; Eisenstein & Loeb 19951 
. Bond & Mvers 1996; Audit et al. 1997; Del Po polo et all 
boOU TShen et al. 2006). Bond & Myers ( 1996|) studied an el- 
lipsoidal collapse model that approximates peaks in the linear 
density field as ellipsoids and accounts for the effects of tides 
on the evolution of overdense patches (see Appendix |B] for a 
more detailed review of the model). Compared with the spher- 
ical case, the key feature of the ellipsoidal collapse models is 
that for a fixed overdensity the collapse epoch will depend on 
the local ellipticity and prolaticity of the shear (or possibly 
density) field. The net effect is that less spherical peaks have 
to overcome additional tidal stretching and require a higher 
overdensity to coll apse. 

SMTOl used the lBond & Myersl(ll996i ) ellipsoidal collapse 
model to derive the dependence of collapse overdensity on the 
ellipticity and prolaticity of the shear field. They found that 
the barrier shape could be approximated by the solution of the 
implicit equation 



^ = l + /3 



(17) 



where (3 and 7 are numerical parameters that must be fit to the 
results of the dynamical model and the squares of the elliptic- 
ity, e, and prolaticity, p, are summed (differenced) if /:> < 
(P > 0). 

Determining the shape of the ellipsoidal collapse barrier re- 
quires a further model for how the typical ellipticity or pro- 
laticity scales with galaxy mass because the barrier shape 
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de pends explicitly on b oth e and p. SMTOl use the results 
of lDoroshkevichI (1 19701) to arrive at a probability distribution 
function for e and p, from which they find the most probable 
prolaticity and ellipticity are /?,np = and e,„p = a j8\fi, re- 
spectively. They then set (5 = 
into Equation ( fTTI ) to find 



(5ec and substituted e^p and 



Be, 



(5ec = (5f 



1 + /3 



cr2(M) 



(18) 



SMTOl found the parameter values /3 w 0.47 and 7 w 0.615 
(m ore recently, values of /3 ~ 0.412 and 7 w 0.618 were found 
by lDesiacquesll2008l) . We refer to the threshold in Equa- 
tion (fTSl l as the "ellipsoidal collapse barrier." For reference. 
Figure [T] compares the ellipsoidal collapse barrier with the 
constant spherical collapse barrier 

SMTOl also suggested an analytical form to approximate 
the first-crossing distribution resulting from the banier (5ec, 
which they calculated numerically using Monte Carlo real- 
izations of the Langevin equation [Equation (O]. They pre- 
scribed the formula 



2 \ '/^ / 2 
v,U{v,) = 2A (1 + v-^'i) ( ^ j exp ( 



(19) 



where q = 0.3, and the constant A = 0.3222 is determined by 
requiring / /ec(i^f)di/,. = 1. 

2.4. The Sheth et al. Barrier and Mass Function 

To improve the agreement between the excursion set mass 
function determined by using Equation (fTsTl and the abun- 
dance of halos in the GIF simulations ( Kauff mann et al.l 19991) 
measured by ST99, SMTOl introduced another parameter, 
asMT ~ 0.707, to modify the ellipsoidal collapse barrier as 



(20) 



and changed the values of the other parameters to /3 = 0.5 and 
7 = 0.6. We refer to Equation ( |20] l as the "Sheth et al. barrier," 
to contrast it with the ellipsoidal collapse barrier because the 
changes are not based on the dynamical collapse model. 

The modified barrier is lower than the Bgc, with the differ- 
ence increasing with increasing mass. In particular, instead of 
converging to the spherical collapse barrier of S^c Sc at the 
largest masses, as expected from the trend towards spheric- 
ity for the rarest peaks, the modified barrier converges to 
SsMT ~ y/osMrSc ~ 0.84^. The value gsmt = 0.707 was justi- 
fied by SMTOl as accounting for the particular choice of the 
halo mass definition in the GIF simulations (identified by the 
Friends-of-Friends algorithm with a linking length b = 0.2; see 
the discussion in §4.1 of SMTOl and Appendix iDb. For ref- 
erence, the Sheth et al. barrier is plotted in Figure[T]alongside 
the spherical collapse and ellipsoidal collapse barriers. 

SMTOl utilized Equation ( l20l i as the collapse barrier to 
calculate an excursi on set mas s function using Monte Carlo 
methods (see also S heth & Torme n 20Q2), and found the re- 
sulting first-crossing distribution to be well-approximated by 
the analytical formula 



I'cfsMTi'^c) = 2A 



X exp 



1 + {\/aSMTVc 



/ 2\l/2 

V 2tt 



(21) 



with flsMT = 0.707, q = 0.3, A = 0.3222. This formula was 
introduced by ST99 as a fit to the GIF simulations, and 
the corresponding halo mass function [calculated from Equa- 
tion ( fT4l )1 is often referred to as the "Sheth-Tormen" mass 
function. The normalization of the first-crossing distribution 
for the Sheth-Tormen mass function is the same as the spher- 
ical and ellipsoidal collapse first-crossing distributions, with 
/ fsMT(yc)diyc = 1- 

3. TESTING THE COLLAPSE BARRIER WITH COSMOLOGICAL 
SIMULATIONS 

As we have noted above, the two main ingredients of the ex- 
cursion set formalism are the collapse barrier used to decide 
which mass elements form halos and the sampling of the pos- 
sible collapse histories via the random walks of the smoothed 
overdensity field. Each of the elements of the excursion set 
approach can be tested against cosmological numerical sim- 
ulations of structure formation. In particular, the correct col- 
lapse barrier, if indeed one can be defined, may be verified 
or falsified using numerical simulations that include all of the 
complications that the excursion set approach aims to circum- 
vent. In the previous section, we described the collapse barri- 
ers motivated by the spherical and ellipsoidal collapse models, 
as well as the modified barrier of SMTOl. In this section we 
compare these collapse barriers with the linear overdensities 
of regions that actually collapse to form halos in numerical 
simulations. 

To calculate the linear overdensity of regions that later 
collapse to form dark matter halos, we utilize cosmological 
simulations perform ed with the Adaptive Refinement Tree 
code (ART Kravtsov "etan |1997|) . One 512^ particle sim- 
ulation models the formation of structure in the Wilkinson 
Microwave Anisotro py Probe (WMAP) Ist-year cosmology 
(ISpereel et al.ll2003L fi„, = 0.3, = 0.7, Qb = 0.04, « = 1, 
(Tg = 0.9, h = 0.7) in a cubic volume L = 250 /!"'Mpc on a side 
(hereafter, the L250 box). This simulation has a gravitational 
force resolution at the highest level of refinement of eL250 = 
7.6 /!"'kpc, and a particle mass of my,.L250 = 9.69 x 10^ Ii'^Mq. 
An additional 10 24^ particle simula tion of the WMAP 3rd- 
year cosmology jSpergel et al.l 12007 . Vl,„ = 0.27, Q\ = 0.73, 
fib = 0.047, n = 0.95, erg = 0.79, h = 0.7) with a larger volume 
(L = 1000 /!"'Mpc; LIOOOW) is used to probe rare objects. 
This larger-volume simulation has a comparably coarser res- 
olution (eiooow = 30 /z-'kpc, »v,Liooow = 6.98 x 10>° r^Mg) 
than the smaller simulation. 

The simulations analyzed in this work were recently used 
bv lTinker et alj (2008) as part of their study of the universal- 
ity of the halo mass function, and we use their halo catalogues 
when calculating th e properties of the dark matter halo popu- 
lation. [Tinker et alj (I2008i) identified dark matter halos using a 
modi fied spherical overdensity algorithm (e.g., Lacey & Co3 
1 19941) . as detailed in their §2.2. Halo membership at a given 
redshift was determined by identifying peaks in the density 
field and assigning dark matter particles to peaks until the 
maximum radius Ra of each halo contains a mean physical 
density of A = M/[4Trp,„(z)R\/3], where M is the sum of the 
particle masses. Results for other halo mass definitions are 
examined in AppendixiDl 

Regions in the linear density field that collapse to form ha- 
los are selected by identifying particles from the z = halo 
catalogue in the simulation v olume at ear l y tinie s. For each 
simulation volume, we use the lZerdovichl(U97(l approxima- 
tion to re-scale the density field at the initial epoch of the sim- 
ulation (z > 50) to a sufficiently early epoch at which the den- 
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8. 150 




50 100 150 200 250 
y [h"' Mpc] 



50 100 150 200 250 
y [h"' Mpc] 



Fig. 2. — Linear overdensity field as a function of smootliing scale for the L250 simulation at a = 0.01. The upper left panel shows the linear overdensity 
field interpolated onto a 512' grid using the cloud-in-cell method (i.e., at full resolution of the simulation); the effective smoothing of this grid corresponds to 
ct(M) ~ 4. The other three panels show the field <5(x,Rw) smoothed on larger scales Rw using a real-space tophat window: Rw = 1.4/j^'Mpc [a(M) = 2.38, M = 
1.0 X 10'2r'],i?w = 2.9r'Mpc [o-(M) = 1.69,M = 8.2 X lO'^r' « M*], and Rw = 6.7r'Mpc [o-(M) = 1.02,M = 1.0 x lO'^r']- 



sity field can be safely considered to be linear We choose 
to re-scale all simulations to the scale factor a = 0.01. The 
density field for each simulation is calculated by a cloud-in- 
cell interpolation of the particle distribution onto a 512^ grid. 
This interpolated density field has a mass resolution equal to 
the particle mass for the L250 box, and eight times coarser 
mass resolution than the LIOOOW box. Hence, the field is 
effectively smoothed on a scale of amaxiM) « 4.2 for L250 
and CTmaxCA/) w 2.1 for LIOOOW. We will use each box only 
on scales a(M) < a,mxiM). Note that since we are interested 
in the regions that form collapsed halos in the Tinker et aD 
(|2008) catalogue, the largest a{M) of interest corresponds to 
the smallest halo for each simulation and all such regions are 

a(M) < (Tn,ax(M). 

3.1. Smoothing the Linear Overdensity Field 

The linear overdensity field smoothed on a scale /?w can be 
calculated directly from Equation ([U by convolving the den- 
sity field with the window function W(x,7?w)- A much more 
computationally-efficient approach is to perform the convolu- 
tion via multiplication in Fourier space to obtain the transform 
of the smoothed overdensity field as 

5R^{k) = kk)W(k,R^), (22) 
and then perform the inverse Fourier transform to arrive at the 
smoothed overdensity ^fi^^,(x). Here, 5 is the Fourier transform 



of the unsmoothed density field and W{k,R^) is the transform 
of the window function [for instance, W{k,R^) = W, for a real- 
space tophat, see Equation (O]. For each simulation volume, 
we compute this convolution for 150 smoothing scales from 
w L/10 to /?w ~ L/256. This results in 512^^ w 134 mil- 
lion overdensity trajectories with 150 steps in a(M), as the 
change of overdensity with decreasing smoothing scale at the 
location of each of the grid cells is equivalent to integrating 
Equation Q with correlated large-scale modes. We have ver- 
ified that the root-mean-squared overdensity fluctuations in 
each box are ((5|^^(x)) = (t^(/?w) simply by averaging over the 
gridded density field, in concordance with Equation (|2]l. 

Figure 12] shows the linear overdensity field 5{a = 1) = 5{a = 
0.01)D(fl = 1)/D(fl = 0.01) of a thin sHce through the L250 
simulation volume. Shown are the unsmoothed field (upper 
left panel) and the field smoothed on scales of a{M) = 2.38 
(/?w = 1.4/i"'Mpc, M w 1O'^/z"'M0, approximately the mass 
of the MiUcy Way halo, upper right panel), a{M) = 1.69 
{R^ = 2.9/i"'Mpc, M = = 8.2 X 1O'^/i"'M0, the present 
collapse mass scale, lower left panel), and <t{M) = 1 .02 (/?w = 
6.7/!"'Mpc, M w 1O'"^/i~'M0, the mass of a large group, lower 
right panel) using a real-space tophat filter [Equation d?)]. 
The figure illustrates a variety of properties of the overden- 
sity distribution and the filter function as the smoothing scale 
is varied. The real-space tophat filter is broad in Fourier 
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Fig. 3. — Smoothed linear overdensity 5. extrapolated to z = 0, as a function of smoothing scale cr(M) for regions that collapse to form halos by z = 0. The 
circles coiTespond to the mean overdensities and the diamonds show the median overdensities, while the errorbars indicate the halo-to-halo scatter. The en'or 
on the mean is significantly smaller than the scatter in all cases. Shown for comparison are the spherical collapse barrier (5c, blue dashed line), the S heth et alj 
1200 Ih ellipsoidal collapse barrier ((5cc, red dashed line), and the collapse barrier associated with the Sheth & Tonnen 1 1999) mass function (green dashed line). 
The upper panel shows the variation in the scatter of barrier heights, where he have used the Greek letter "S" to denote this scatter [not to be confused with 

<T{M)l 



space, so the variations in overdensities across many intervals 
in the smoothing scale are correlated. A sharp A:-space filter 
would tend to decorrelate the overdensities on small smooth- 
ing scales from larger scales as independent frequency modes 
are added with increasing a(M). The largest fluctuations in 
the unsmoothed overdensity field are identifiable across the 
smoothed fields, reflecting the relation between the initial 
density fluctuations and the eventual formation of massive 
structures. 

3.2. Smoothed Overdensities of Collapsed Regions at z = 

We can use the smoothed density fields to connect the fi- 
nal mass of a collapsed region [and its associated fluctuation 
scale (t{M)] with its initial smoothed overdensity linearly ex- 
trapolated to the epoch of observation. For each halo iden- 
tified in our catalgoues at z = 0, we calculate the center-of- 
mass of the halo particles from their positions in the lin- 
ear density field at a = 0.01 and use the window-smoothed 
field to compute the overdensity within the lagrangian radius^ 

* Although this is not entirely self-consistent (the smoothing scale is re- 
lated to the mass scale by M oc (1 +(5)R^), it is clear that the error is small if 
this calculation is done at an epoch when S is very small. We have tested that 



R= [3M/47rp„,]'/"' about this location. This overdensity is 
then linearly extrapolated to z = to serve as an estimate of 
(5m(x). We have checked that all of our conclusions are robust 
to specific choices regarding the smoothing procedure, such 
as the choice of initial positions (the x) of halos in the initial 
density field and the range of smoothing scales. 

Figure |3] shows the distribution of such smoothed linear 
overdensities extrapolated to z = as a function of a{M). The 
shaded regions represent the probability distribution p{5,a) 
for regions that collapse to form halos defined relative to 
a A = 200 spherical overdensity identified in the LIOOOW 
(601,448 halos) and L250 (73,720 halos) boxes. The median 
(colored diamonds) and mean (colored circles) of 5{z = 0) in 
bins of width A(7(M) = 0.25 are measured for the halo pop- 
ulation and shown for comparison. The distribution at fixed 
<t{M), p(5\<t), is approximately log-normal in shape, with a 
width that scales as I]((5) w 0.3ct' *'(M) (indicated by the error 
bars in Fig. O. The scatter in 6{z = 0) at fixed (j{M) reflects 
both the intrinsic scatter in the linear overdensity of collapsed 
regions and the limitations of our method to measure 6{z = 0) 
reliably for any individual halo. The error on the mean or 

the epoch we use a = 0.01 is sufficiently early for this purpose. 
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Fig. 4. — First-crossing distributions for c ommon collapse barri ers from tlie literature. The upper row shows distributions calculated with th e sharp i-space 
tophat filter for the spherical collapse barrier IGunn & G otf 1972". gray line, upper left panel), the ellipsoidal collapse barrier ll Sheth et"aI]|2001L red line, upper 
middle panel), and the collapse barrier associated with the Sheth & Tormen 1 1999) mass function I Sheth et al. 2001, green line, upper right panel). The colored 
lines indicate the sharp ir-space excursion set first-crossing distributions calculated using the Zhang & Hui 1 2006) analytical method. The histograms indicate 
Monte Carlo realizations of the first-crossing distributions calculated by integrating the Langevin Equation (5). The dashed black lines indicate the analytical 
form for the Press & Schechtei ( 1974) mass function (upper left panel), the S heth et ah 12001,) fit to the ellipsoidal collapse first-crossing distribution (upper 
middle panel), and the Sheth & Tormen 1 1999) mass function (upper right panel). Also sho wn is t he Tinker et al. 1 2008) mass function determined from a large 
suite of cosmological simulations (blue line, all panels). In all cases, the Monte Carlo and [Zhang & Hui 12006) first-crossing distributions agree well. Since 
the ellipsoidal collapse banier hes above the spherical collapse barrier (see Fig ure [T}, th e corresponding first-crossing distribution predicts fewer galaxy mass 
{uc > 0.5) halos than does the Press & Schechter, l_1974) formula. The Sheth & TormeiJ (1999) mass function produces more high-mass halos by lowering the 
associated collapse barrier below the spherical collapse barrier Bsc = <5c at small a{M} (see Figure [T). For comparison, the bottom row shows the first-crossing 
distributions for the same barriers calculated using a real-space tophat filter (histograms, bottom row). In each case, using a real-space tophat filter produces 
fewer halos in the excursion set calculation than does the sharp i:-space tophat for most halo masses. 



median in any aiM)-bm is much smaller than E((5). 

At all measured halo masses, the mean and median lin- 
ear overdensity of the halo population at 5{z = 0) exceeds 
the spherical collapse overdensity 5c. The dependence of the 
mean and median overdensity on a{M) measured in simula- 
tions increases in a manner that resembles the functional form 
of the ellipsoidal collapse barrier (5ec [Equation (fTsl )] and mod- 
ified Sheth et al. barrier (5smt, but with a different normal- 
ization. Importantly, at the lowest values of variance (largest 
masses) probed by the simulations {a(M) ^ 0.5] the measured 
overdensities of the collapsed objects are larger than both the 
modified Sheth et al. barrier (5smt cind the spherical collapse 
barrier 5c. We note that while our L = l/i"'Gpc ACDM sim- 
ulations do not probe the highest masses (M > 5 x lO'^M©) 
and the the largest scales (cr(M) < 0.4) imaginable, in Ap- 
pendix|D]we show that the conclusions developed from these 
results do not change even if the characteristic overdensity of 
regions that collapse to form halos does not asymptote exactly 
to 5c for cr ^ 0. In the next section of the paper, we explore 
the implication of these results for halo mass functions calcu- 
lated using the excursion set formalism. 



4. EXCURSION SET FIRST-CROSSING DISTRIBUTIONS AND MASS 
FUNCTIONS IN COSMOLOGICAL SIMULATIONS 

The results presented in the previous section indicate that 
the linearly-extrapolated overdensities of regions that collapse 
in cosmological simulations behave in a manner analogous to 
the expectations of the ellipsoidal collapse model. In this sec- 
tion, we test the second ingredient of the excursion set ansatz: 
the calculation of the first-crossing distribution and associated 
h alo mass function w ith a particular barrier 

[Tinker et al.l (l2008l) used a large suite of cosmological simu- 
lations to determine an accurate numerical fit to the abundance 
of dark matter halos as a function of their mass. They found 
that the first-crossing distribution that corresponds to the halo 
mass function measured in simulations can be well-described 
by the function 



e-YVc 



exp 



" 5} 



(23) 



For A = 200 spherical overdensity halos, the best fit mas func- 
tion parameters are Aj = 0.482, dj = 1.97, ej = 1, ,?t = 0.51, 
and h j = 1 .228 (with x^/j^ =1.14, see table C4 of lTinker etal] 
l2008h . 
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Fig. 5. — Renormalized first-crossing distributions for ellipso idal collapse 
vs. th e [Tinker et aTl 120081) simulation results. The form of the [Tinker et alj 
120081) mass function is normalized such that / fj(uc)duc = 1, but the 
simulation results only probe masses Uc ^ 0.5 that contain roughly 60% 
the mass of the universe. This figure demonstrates that any associated 
uncertainty with the normalization of the Tinker et alj t2008i) mass function 
(solid line) cannot reconcile the difference with the elhpsoidal coll apse 
first-crossing distribution. If the ellipsoidal collapse barrier (Equation 1 18) 
first-crossing distribution is renormalized to match the Tinker et al. (2008) 
mass function at Uc = 3 (dashed line), or renormalized to match the integral 
of the Tinker et al. 12008) mass function alUc > S (dotted li ne), the resulting 
distribution is still discrepant from the [Tinker et al.l j2008[) results. Hence, 
the shape and normalization of the ellipsoidal collapse barrier excursion set 
mass function differs from the halo mass function determined by N-body 
cosmological simulations. 



Figure |4] compares the first-crossing distribution given by 
Equation (|23T l with first-crossing distributions calculated from 
the excursion set formalism with the collapse barriers de- 
scribed in § |2] Where possible (i.e., in the case of sharp k- 
space tophat window function) we have checked our Monte 
Carlo calculation of the first-crossing distribution against the 
direct solution of th e Volterra equation using the method of 
IZhang & Hull (l2006i see Appendix and have found excel- 
lent agreement. 

The spherical collapse model (left column), corresponding 
to the Press & Schechter (1974) mass function, displays the 
well-known deficit of massive halos and overabundance of 
low-mass halos compared with simulations (see the discus- 
sion in § [2.2[ and Appendix |A|l. As demonstrated first by 
iBondet all (119911 see their Figure 5), changing the window 
function from the sharp ^-filter to a real-space tophat filter acts 
to reduce the abundance of galaxy-mass halos rather than im- 
prove the agreement between the spherical collapse and sim- 
ulation first-crossing distributions. 

Using the ellipsoidal collapse barrier presented by SMTOl 
(see Equation [TsT l with the excursion set calculation results 
in a lower first-crossing distribution at galaxy masses (Fig- 
ure |4] middle column) compared to the spherical collapse 
model. These first-crossing distributions lie slightly below the 
SMTOl fitti ng function [Equation (fT9]l1, and well below the 
[Tinker et"ai . (2008) simulation results for > 1- The ellip- 
soidal collapse barrier converges to the spherical collapse bar- 
rier as cr(M) 0, but at the mass scales probed in simulations 
[o-(M) > 0.5, M < 1O'^/!"'M0] the ellipsoidal collapse barrier 



is considerably larger (6^^ }^ 1-8). The larger barrier height 
tends to suppress the abundance of halos in this mass range 
as calculated by the excursion set formalism . Similar results 
have been obtained by other authors (e.g., San dvik et al.l200'7l 
their Figure 3), but the convergence to the spherical collapse 
model at i^c >> 1 is seldom commented upon. For consis- 
tency between the treatments of spherical and ellipsoidal col- 
lapse, the abundance of halos calculated for small ellipticities 
and prolaticities (i.e., at large masses) must be the same in 
both models. As with the spherical collapse model, the ellip- 
soidal collapse excursion set mass function calculated with a 
real-space tophat filter acts to lower the abundance of galaxy- 
mass halos (Figure,!?] bottom, middle panel). 

The Sheth et al. baiTier [Equation d^Ull], which SMTOl pre- 
sented as the barrier corresponding to the GIF simulation 
(ST99) mass fun ction, produc es a first-crossing distribution 
that lies below the lTinker et al.l (2008) simulation results (Fig- 
ure |4] right column). Note also that the predicted distribu- 
tion does not agree with the SMTOl first-crossing distribu- 
tion (Equation (l2Tl i. dashed line. Figured upper right panel), 
which this barrier is d esigned to describe . The relative dis- 
agreement between the lTinker et al.1 (l2008h and SMTOl first- 
crossing distributions should be noted. This disagreement, al- 
though relatively small on the scale o f this figure, is signi ficant 
in terms of the halo abundance (see [Tinker et al.ll2008b . For 
the first-crossing distribution determined using a real-space 
tophat window function (Figure |4] lower right panel), the 
Monte Carlo calculation lies below the sharp A:-space window 
function results and is discrepant with the simulation result. 

Since all simulated mass functions are limited by their par- 
ticle resolution one might wonder if the differences with the 
simulations owe to the differences in normalization, as all of 
the excursion set first-crossing distributi ons are normalized 
such that / f{vc)dvc = 1 . The form of the [Tinker et al.1 ([20081) 
first-crossing distribution [Equation ( i23T ll that is plotted in 
FigurelHis constructed to require that j fj{vc)dvc = 1. How- 
ever, the free parameters of this function are fitted in a regime 
that incorporates only « 60% of the available mass. To ad- 
dress this question. Figure |5] shows the ellipsoidal collapse 
first-crossing distribution normalized at a fixed mas s scale 
(instead of the integral constraint) to match the Tinke r et"an 
(2008) first-crossing distribution at large masses {Vc = 3), or 
constrained to provide the same integral over the same mass 
range (at Vc > 3). The functions clearly differ in shape and 
not simply in normalization. 

Of the three barriers we examined, the excursion set first- 
crossing distributions calculat ed from the Sheth e t al. barrier 
most closely approximate the Tinker et al.l (l2008l) simulation 
results. This result is not surprising, given that the ST99 func- 
tion is itself a fit to cosmological simulations. The ellipsoidal 
collapse barrier produces a first-crossing distribution that dif- 
fers substantially from the ST99 fitting function, even though 
this mass function is frequently associated with ellipsoidal 
collapse model. 

The results presented in this section demonstrate that ex- 
cursion set predictions disagree with ACDM simulation mass 
function results for all of these barriers. Given that the over- 
densities of the collapsed regions behave similarly to the bar- 
rier shape expected for the ellipsoidal collapse, our results im- 
ply a manifest failure of the excursion set ansatz as a method 
of computing the abundance of collapsed objects. In Ap- 
pendixini we show that our results hold for other halo defini- 
tions, including Friends-of -Friends halos and other spherical 
overdensity definitions. 
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5. DISCUSSION 

In the excursion set theory, the abundance, formation time, 
and bias of dark matter halos are assumed to be directly linked 
to the initial linear overdensity field extrapolated to a given 
epoch using linear growth rates. This assumption significantly 
simplifies modeling, but clearly needs to be tested against di- 
rect numerical simulations. Our results show that the linear 
overdensities around the Lagrangian positions of the centers- 
of-mass of collapsed halos in simulations behave in a man- 
ner analogous to the collapse overdensity barrier predicted by 
SMTOl using an ellipsoidal collapse model. The similar be- 
havior suggests that this model captures the main physics be- 
hind the nonlinear gravitational collapse around peaks in the 
initial density field. 

At the same time, the failure of the excursion set ansatz to 
predict correct abundance of collapsed halos with the ellip- 
soidal collapse barrier demonstrates that the ansatz is flawed. 
Several aspects of the ansatz may be responsible for its fail- 
ure to accurately describe simulation results. For example, 
assumption that each mass element can be assigned to a col- 
lapsed halo using only the local overdensity independently 
of its environment is definitely problematic. The threshold 
overdensity for collapse (the barrier) in either the spherical 
or ellipsoidal collapse models is predicted for the volumes 
centered on a density peak, and not for a random mass ele- 
ment in the field. While the actual collapse will occur around 
the density peaks, additional mass near a peak may collapse 
onto it even though the extended region may not satisfy the lo- 
cal collapse condition. Our study highlights the failure of the 
excursion set ansatz on a statistical basis for the entire halo 
population, but a physical model for the collapse and forma- 
tion of da rk matter halos m ust also succeed object-by-object 
(see, e.g.. iKatz et al.|[T993h . The methods for calculating the 
abun dance of halos that treat individual peaks i n the density 
field (iBond & Myerslll996l: iMonaco et al.ll2002l Dalai et al., 
in prep.) may therefore afford a better way to predict the for- 
mation, masses, and abundances of collapsed objects. 

The shape of the effective collapse barrier is determined 
both by the distribution of overdensities in a smoothed Gaus- 
sian field (e.g., with ellipticities and prolaticities, or other 
properties, sufficient for collapse) and by the dependence of 
the collapse condition on the overdensity (and/or other prop- 
erties). The halo formation time can be defined relative to 
the mean assembly history (e.g., Wechsler et al. 2002), but 
the specific choice of definition can influence the relation be- 
twe en formation t ime and halo mass (for a recent discussion, 
see iLi et ani2008h . In terms of the linear overdensity and the 
effective barrier B, one definition of halo collapse is simply 
that of Eq. (fT2l i. Of course, determining the effective barrier 
B or measuring the collapse epoch Zc complicates the matter. 
The definition of collapse can be connected to the physical 
properties of halos through dynamical models, but identify- 
ing evolutionary phases in such simple models with the ac- 
tual nonlinear growth of dark matter halos may be incorrect 
or inaccurate. For ins tance, the ellipsoidal collapse model of 
iBond & MyersI (11996 ) associates halo formation with the col- 
lapse of the longest ellipsoid axis, while freezing the collapse 
of the shorter axis at a particular point to prevent a density sin- 
gularity. Halo virialization, however, may be associated with 
the collapse of another axis or condit ional pro perties of shape 
of the density or s hear fields (e.g.. iLee & Sh andarini ll998t 
[Monaco et al.ll2002h . These issues warrant a more careful ex- 
amination that we defer for future work now that we have a 



statement of the problems at hand. 

Dark matter halos reside in special locations of the density 
field, and have d ifferent cluste ring properties than bulk mat- 
te^Xe^., Hiiii 1984; Efstat hiou et al.lll988l: ICole & Raised 
'1989b). Models of halo bias are tightly conn ected with the 
effective collapse barrier. iMo & Whitd (1 1996) presented the 
idea that the biasing of halos is connected to the rate at which 
density trajectories cross two separate barriers, and used this 
concept to calculate the Lagrangian bias of halos in the Press- 
Schechter model. SMTOl adapted these ideas to calculate the 
halo bias implied by their modified ellipsoidal collapse model. 

In principle, the same collapse barrier should predict con- 
sistently both the bias and mass function of halos. How- 
ever, while the ST99 mass function models the halo abun- 
dance reasonably well, ISeliak & WarrenI (|2004|) found that 
the SMTOl bias model does n ot work well a t low m asses. 
For mas sive halos at high V c, ICohn & Whitd (l2008l) found 
that the iMo & Wh"it3 d 19961) bia s scaUng works better than 
the SMTOl model (however see iHu & Kravtsovll2003h . In- 
terestingly, 'Ree d et aH (l2008l) found that early-forming halos 
in their simulations were well-modeled by SMTOl at peak 
heights as large as z/^ ~ 4 (at small Vc, the large-scale bias they 
measure drops b elow the SMTOl model , in a manner similar 
to that found bv ISeliak & WarrenI l2004h . We note that if the 
effective collapse barrier converges to the spherical collapse 
overdensity 5c for large mass or highly-bi ased halos, then one 
might expect the halo bias to mimic the lMo& Whitd (11 996h 
scaling at large Vc- We speculate that the intriguing devia- 
tions between the models for halo bias and abundance of halos 
are connected to the current discrepancy between the effective 
collapse barrier and the abundance of halos in the excursion 
set ansatz that we discuss in this paper. 

Lastly, the effective collapse barrier may be connected 
with the assembly bias phenomenon wher e dark matter halo 
clustering correlates with formation time. iGao et all (l2005h 
found that low-mass halos (M < M^) that formed early were 
much more strongly cl ustered than late-forming halos of 
the sa m e mass (see also Sheth & Tormenll2004l : iHarker et alj 
l2006h . IWechsler et aP (l2006h demonstrated that the sense 
of assembly bias reverses at high-mass (M > M^) such that 
late-forming halos were more strongly clustered than early- 
forming halos at fixed mass, and that these correlations are 
reflected in the relative bias of halos with different concen- 
tratio ns at fixed mass (see also IWetzel et alJl2007t iJing et all 
|2007|) . 

Zen tneij ( |2007|) in his review showed that window functions 
that are local in real-space rather than Fourier space naturally 
result in early-forming halos that reside in underdense re- 
gions, as reported for high-mass halos (M > M*) in the simu- 
lations. The standard implementation of the excursion set for- 
malism assumes that the process of nonlinear collapse can be 
encapsulated into the assignment of the collapse barrier. For 
high-mass halos, tidal influences and nonlinear interactions 
with nearby objects should be minimal because halos with 
masses significantly larger than form from nearly spheri- 
cal peaks and usually dominate their local environments. As 
a result, the excursion set assumption should be most valid 
for high-mass (M ^ M*) halos and leads to a natural picture 
where early-forming, high-mass halos become less strongly 
clustered than their late-forming counterparts. 

The reversal of the environment-dependent halo forma- 
tion trend at low-mass may owe to the truncati on of small 
halo g r owth by nearby str ucture as suggested bv i Wanp et alj 
( I2007t) . iDalal et alj ( 120081) greatly extended this work, vaUdat- 
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ing the high-mass trend in a set of scale-free numerical sim- 
ulations and showing that environmental influences on halo 
bias, concentration, and formation time at fixed mass could be 
accounted for by considering the "peak curvature", dSR^/da, 
in addition to peak height. The peak curvature serves as some 
proxy for environment as peaks with greater curvature lie in 
relatively underdense environments. iDalal et aTl (l2008l) also 
used a toy model to demonstrate that at s mall masses environ - 
mental effects, such as those suggested bv lWang et al.l (120071) . 
can truncate halo growth and drive early-forming halos to be- 
come less anti-biased as they are advected by the larger-scale 
matter field. 

Our results relate to these findings by providing a new out- 
look on the connection between halo abundance, the effective 
collapse barrier, and the excursion set formalism. If the ex- 
cursion set formalism fails to account properly for the abun- 
dance of halos using an appropriate form for the collapse 
barrier, then it may also fail to describe reliably the con- 
nection between formation time, mass, and properties of the 
density field. The effective collapse barrier may well be a 
funct ion of additional parameters beyond the local density 
(e.g., IChiueh&Lejl200H ISandvik et al. l IIoOTI) . may incor- 
porate informationabOT (e.g.. lZentneii 
I2007t iDesiacquesI [20081: iDalal et allllOOSh . and may require 
additional parameter dependencies that to account for the 
nonlinear collap se of overdensities (e.g., Wang et al. 2007; 
IDalal et al.ll20^ . 

At present, the excursion set theory provides the main 
framework used to develop heuristic understanding of simula- 
tion results and to formulate fits to simulation results for halo 
abundance, clustering, and other halo properties. The method 
succeeds in a gross sense. Excursion set theory identifies the 
fundamental scale in the problem, the mass where a(M) ^ S^.. 
Below this characteristic mass the halo abundance per unit 
mass has a simple, power-law form, while above this mass the 
halo abundance drops rapidly. However, a precise understand- 
ing of halo abundance and clustering beyond the gross accu- 
racy provided by the excursion set ansatz is now necessary. 
Contemporary and forthcoming efforts to use measurements 
of the abundance and clustering of galaxy clusters to constrain 
cosmological parameters, as well as comprehensive statistical 
studies of galaxy formation and evolution, only highlight the 
need for a sound understanding of halo abundance and assem- 
bly, both globally and as a function of environment. We need 
to explore amendments and alternatives to the excursion set 
model to make progress in our understanding of halo forma- 
tion. 

6. SUMMARY 

In this paper we have presented tests of the excursion set 
ansatz against cosmological simul ations. Usin g a subset of 
cosmological simulations from the [Tinker et al.l (I2008I) study 



of the halo mass function, we identify the locations in the lin- 
ear overdensity field that later collapse to form dark matter 
halos. We demonstrate that the dependence of the linear over- 
density of these regions on mass or smoothing scale a{M) re- 
sembles predictions of the ellipsoidal collapse model. While 
the effective collapse barrier of simulated halos behaves anal- 
ogously to the simple ellipsoidal collapse barrier, the simu- 
lated halo mass function is inconsistent with what the excur- 
sion set ansatz predicts for such a barrier This inconsistency 
implies that the excursion set ansatz is not valid and cannot 
be used reliably to predict halo abun dance or bias . 

The modified collapse barrier of ISheth et aP (1200 Ih dif- 
fers significantly from the physical behavior calculated by 
the ellipsoidal collapse model, which for example predicts 
convergence to the spherical barrier (5 6c) for the rarest 
peaks (cr 0). In vi ew of this, the interpretation of the 
ISheth & Tormeiil ([19991) mass function as a prediction of the 
ellipsoidal collapse model for the abundance of dark matter 
halos is not correct. 

The impressive statistics of ever-larger dark matter simula- 
tions will likely continue to uncover increasingly subtle varia- 
tions on the classical picture of dark matter halo formation. 
In this work, we identify a striking inconsistency between 
the effective collapse barrier of simulated halos and excur- 
sion set formalism predictions for their abundance. Our re- 
sults also demonstrate that there is still much to learn and 
understand about the conditions for and the process of halo 
collapse, which warrants further studies that critically revisit 
these issues using modern large cosmological simulations. 
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APPENDIX 
A. SPHERICAL COLLAPSE 

A simple approximation for the dynamical evolution of an overdense region is the spherical collapse model. The dynamical 
equation for the expansion an d collapse of a spher ical overdensity A of physical size ^phys(z) in a flat universe with a cosmological 
constant can be written (e.g. JGunn & Got3ll972h 

= |7rGpAflR-^7rGpmaR[l + A(0], (Al) 
d/^ 3 3 

where or = Rphys/R is the scale factor of the region, pA = ^aPc is the dark energy density, and J^a is the dark energy density 
parameter. The corresponding growing linear overdensity 6 obeys the differential equation 

d^S dS 

-—r + 2H—= AttGp„5 (A2) 
df ^ at 

(e.g.. lLifshitzlll946t lPeebleslll965h . while the physical overdensity evolves as 

A(f)=^-1, (A3) 
where a is the universal scale factor. The initial conditions for evolving Equation |AT] are simply 

flR(rinit) = fl(finit)(l - '5(rinit)/3), (A4) 
flR(finit) = //(finit)flR(finit)-a(finit)//D(finit)(5/3, (A5) 

where Hq = D/D describes the rate of change of overdensities. Initially A(f) ~ 5{t), as can be checked by Taylor-expanding 
Equation IA3 1 and comparing it with Equation I A4I but eventually the quantities diverge as the overdensity begins to exceed the 
applicability of the linear order approximation. As the spherical region begins to break from the universal expansion and reaches 
a maximum radius at the turn-around time f = fta, the physical overdensity reaches 1 -I- A(/ = fta) — 5.55 while the linear overdensity 
is 6{t = fta) ~ 1 .06. As the region collapses to a point of zero size, A ^ cx3 while the linearly extrapolated overdensity approaches 
a value of 

Bsc = '5sc = 5(fmit)£>(Zc)/£>(zinit) = '5c (Spherical Collapsc), (A6) 

where 5c = 1 . 686 is often called the linear collapse overdensity for the growth of spherical perturbations (Gunn & Gott"1972t 
lPeebles|[T980l) for an n,„ = 1 universe (S^ has a weak dependence on cosmology, see, e.g.. Bond & Mvers 1996; Eke et al. 1996). 
If the spherical collapse of the region to zero size is associated with the formation of a dark matter halo, then Equation |A6] can 
be used as the "spherical collapse barrier" for purposes of calculating the first-crossing distribution associated with spherical 
collapse. 
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B. ELLIPSOIDAL COLLAPSE AND ITS MODIFICATIONS 



The gravitational co llapse calculation has been generalized to model nonspherical collapse by a number of authors (e.g., 
ZelMovich' 4970llNa riai &Fuiimoto 1972; Hoffman 1986; Bertschi nger & Jainlll994t lEisenstein & Loebll 19951; iBond & MversI 
1996: Audit et al.l]r997: Del Popolo et al. .2001: Shen et al.. .2006) . Below we focus on an ellipsoidal collapse model by 



Bond & Myers! (119961 see their §2.1.3 and Appendix A), which treats the gravitational col l apse o f a homogeneous ellipsoid 



by separately following the coupled evolution of each axis a, of the ellipsoid. ' Bond & MversI (1199 6') showed that the dynamical 
equation for the expansion and collapse of an ellipsoidal region with linear overdensity 6(t = fink) = Sq at an initial time finit can be 
written 



8 

= -7rGpAfl/-47rGp,„fl; 



where the index / = 1,2,3 indicates a principal axis. The term 

fc,- = -- + aifl2fl3 / 
-5 Jo 



dr 

(fl2 + r)n,i^,(fl2+r)i/2 



(Bl) 



(B2) 



describes triaxial contributions to the gravitational acceleration. A linear order approximation of the effects of external tides on 
the evolution of the region are included through the factors AJ(0 = A, - (5/3, which are written in terms of the eigenvalues of the 
strain tensor 

A3 = -(l + 3e + p), (B3) 



■2/7), 



■3e + p). 



(B4) 



(B5) 



Here, e > is the ellipticity, -e < p < e is the prolaticity, and the eigenvalues are ordered A3 > A2 > Ap This linear order 
approximation to the effects of tides grows as A,(f) cx 6(t) cx D(t). The physical overdensity in the ellipsoidal collapse model 
evolves simply as 

A{t)=— 1, (B6) 



where a is the universal scale factor The initial conditions for evolving Equation IB II are 

Oiihmd = ^2(fi„it)[l — A,(finit)], 
diitinit) = ^^(finit)fl/(finit) ~ <2(finit)//£)(finit)A,(finit), 



(B7) 
(B8) 



The evolution of the region is then determined by the cosmology, the initial overdensity (5o, the ellipticity e, the prolaticity p, and 
the initial universal scale factor a(t = fi„it). For a spherical system (e = 0, p = 0, oi = 02 = 03), the dynamical equations reduce to 
Equation ( lAlb . 

Sheth et al. (2001) used the results of iDoroshkevichI d 1 9701) to show that, in the context of the lBond & MversI d 19961) elhpsoidal 
collapse model, the formula 



g(e,p\S) = 



exp 



5f 

'2(7 



(B9) 



provides the expected distribution of ellipticities and prolaticities for the shear field of Gaussian random overdensities and corre- 
sponds to a distribution of effective collapse barriers for halos (as a function of e and p) that describes the overdensity at which 
the last (longest) principal axis collapses. Sheth et al. (2001) provided an empirically-determined, imphcit functional form to 
approximate the shape of the ellipsoidal collapse barrier in terms of e and p, 



Sec(e,p) 



= l+(3 



5(e^±p^) 



(52 



(BIO) 



which follows the most probable {p = 0) trend of the collapse barrier distribution well for the parameter values /3 = 0.47 and 
7 = 0.615 (a more recent calculation has found (3 = 0.412 and 7 = 0.618, Desiacques 2008). For the most probable prolaticity 
p = 0, the maximum of the probability distribution gie^p = 0\S) follows the ridgeline Cmp = (ct/S)V5. Substituting p = and 



: Cmp into Equation IB 1 01 yields a characteristic ellipsoidal collapse barrier in terms of the overdensity variance as 



(BID 



Note that this barrier reduces to the spherical collapse barrier [B = 5^ 

reference, the ellipsoidal collapse barrier is plotted in Figure[T]with the lSheth et al.l (l2001l) parameters (dashed red hne). 



E quation (flStl for la rge halo masses (small variances). For 
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C. AN INTEGRAL METHOD FOR FIRST-CROSSING DISTRIBUTIONS WITH SHARP AT-SPACE FILTERING 

Typically, the first-crossing distribution is determined with a Monte Carlo approach (effectively by integrating Equation|3]l. For 
the case of sharp ^-space filtering, where smoothed overde nsity executes a Markovian random walk with the smoothing scale, and 
barriers with a suitably weak dependence on the variance, IZhang & Huil (f2006j showed that by properly accounting for the rate 
of first barrier crossings in an ensemble of trajectories, an integral relation for the first-crossing distribution could be expressed 
in terms of the functional form of the barrier B and the probability P(S,S) that a trajectory 6 first crosses the barrier near S. 
Specifically, in terms of the barrier B(S), the first-crossing distribution satisfies 

rS rB(S) 

1= f{S')dS'+ P{5,S)dS (sharp-;t filtering). (CI) 

Jo J-oo 

The first term on the right hand side accounts for trajectories that have crossed at scales larger than S, while the second term 

P{d,S) = PQ{6,S)- f dS'f{S')Po{d-B{S'),S-S') (C2) 
Jo 

subtracts the rate of down-crossings from the probability 



that a trajectory given crosses the barrier near 5'. Equation ICll can be differentiated and combined with Equations IC2HC3I to 
produce a Volterra integral equation of the second kind for the the first-crossing distribution 

f(S) = gi(S)+ f f(S)g2(S,S')dS' (sharp-^ filtering), (C4) 
Jo 



with 



B(S) ^dB 
S dS 



82(S,S'): 



^dB B(S)-B(S') 



Po(B(S),S) (C5) 



Po(B(S)-B(S'),S-S') (C6) 



dS S-S' 

Hence, given a barrier shape B{S), IZhang & Huil (120061) have provided a helpful method for calculating the first-crossing distri- 
bution for a A:-space filter via Equation lC4l . 

D. THE EXCURSION SET AWSArZ AND HALO MASS DEFINITIONS 

The results of ^demonstrate a disconnect between the characteristic linear overdensity of regions that collapse to form dark 
matter halos and the collapse barrier required to reproduce the abundance of those same halos using the excursion set formalism. 
While the results of ^are internally consistent, one might wonder if the failure of the excursion set ansatz was peculiar to the 
A = 200 spherical overdensity halo definition. In this Appendix, we demonstrate that the excursion set ansatz also fails for other 
common halo definitions (specifically, A = 100 spherical overdensity halos, A = 600 spherical overdensity halos, Friends-of- 
Friends halos, and halos defined by spherical regions of size R = 2/?2oo)- Since these halo definitions span the most practical 
definitions found in the literature, the results of this appendix present an exhaustive study of how our results depend on the halo 
and mass definitions. Further, since the largest simulation we study has a l/!"^Gpc^ volume, the linear overdensity of regions 
with very large mass (M200 5 x 10'^/i~'Mq) are not probed by our simulations. Below, using extrapolations of the 5-a{M) 
trend for various halo mass definitions, we demonstrate that even if 5 ^ 5c as a{M) the excursion set mass functions do not 
reproduce the simulated mass function for any halo mass definition we consider 

To repeat the calculations in ^for other halo definitio ns, we must construct additional halo catalogues. For the spherical 
overdensity halo definition, we follow [Tinker et al.l (l2008l) and define halos with an overdensity A relative to the background 
density p„, as the particles within a radius /?a around density peaks. For the Friends-of-Friends halo definition (e.g., Davis et all 
[1985), we adopt the standard linking length ofb = 0.2. For halos defined by spherical regions of size R = 2R/^, we use the A = 200 
catalogue to identify halos and redefine the halo masses by assigning all particles within 27?2oo of the center-of-mass membership 
in the halo. If the radius R = 27?2oo for one halo includes the center-of-mass of a smaller halo, the smaller halo is discarded from 
the catalogue. 

For each halo definition, the mass function d« / dM is determined by constructing a histogram for the halos by binning in mass. 
The first-crossing distributions /(i^e) corresponding to e ach mass fun ction are calculated using Equation ( fT4] i. For each mass bin 
we calculate jack-knife errors, as described in detail bv lTinker et al.l (2008). For two spherical overdensity definitions, A = 200 
and A = 600, we s i mply adopt a mass function of the form of Equation ( |23] | with the best-fit parameters determined in Appendix 
C of lTinker etakl (l200l . For A = 200, At = 0.482, dj = 1.97, ej = 1, gj = 0.51, and hj = 1.228. For A = 600, At = 0.494, 
dj = 2.56, ct = 0.93, ^t = 0.45, and /jt = 1.553. These analytical mass function fits make use of the wide range of simulations 
studied by lTinker et a l. (2008). For mass functions for the other halo definitions, we rely on our halo catalogues for the LIOOOW 
simulation and represent the mass function with binned values and uncertainty estimates constructed from these catalogues. We 
have checked that the first-crossing distribution calculated from the binned mass function for the A = 200 and A = 600 halos in 
the LIOOOW simulation match the corresponding analytical fits from lTinker et al.l (l2008h extremely well, and we therefore expect 
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that the binned mass functions and first-crossing distributions for the other halo mass definitions are reliable estimates of the halo 
abundance over the mass range probed by the LIOOOW simulation. The LIOOOW simulation is sufficient for our needs, as the 
primary constraint comes from the most massive halos in our catalogues. 

We construct the distribution of linear overdensity 5 as a function of the smoothing scale a{M) for the regions that collapse to 
form halos in the manner described in ^3.21 using the same set of smoothed linear density fields calculated in ^3.11 The mean of 
the S{(t) distribution is calculated in three bins of width Act = 0.25. For the A = 200 halo definition, the calculation results in the 
distribution of S vs. a{M) shown in Figure[3](at a{M) < 1.2). 

To calculate an excursion set mass function from the overdensity distribution for each halo mass definition, the mean overden- 
sity (5 as a function of a{M) is fit with two analytical forms to produce two model collapse barriers. We first fit the function 

Snt = Sc[l + f3(cTyS^^r], (Dl) 

used by lSheth et al.l (1200 Ih to represent ellipsoidal collapse, allowing (3 and 7 to vary. By construction, this function converges 
to (5fit — * 6c as (7 ^ 0. We also fit a linear function, 

Si,,=Aa + b, (D2) 

allowing A and b to vary. The intercept b in general is smaller than Sc', this functional form thus allows us to test how results 
would change if the barrier does not asymptote to Sc for low values of cr. The best fit parameters for the effec tive coll apse barrier 
for each halo mass definition are reported in Table[T] We the n use the best fit p arameters for Equations ID 1 1 and ID2 1 to calculate 
a sharp ^-space first-crossing distribution via the method of IZhang & Huil ([2006) [i.e.. Equation ( IC4l i1. and compare with the 
simulated first-crossing distribution provided by the halo catalogue. 

Figure|6a]shows the distribution of linear overdensity 6 with smoothing scale cr(M) for the A = 200 halo catalogue. As in Figure 
[3] the mean overdensity of regions that collapse to form halos in the A = 200 catalogue lie above the spherical collapse barrier 
Sc- Analytical fits to the mean of the overdensity distribution (i.e., the effective collapse barrier for this halo mass definition) 
show that a simple linear extrapolation possibly suggests that S ^ 1 .5 as <t ^ 0. The excursion set first-crossing distributions 
calculated for models of the mean overdensity in regions that form A = 200 halos show that the differences in the fits at cr < 0.45 
have little influence on the resulting halo abundance, as the la rger 6 > 6c at scales cr > 0.45 suppresses the abundance of halos at 
I'c < 4 relative to the spherical collapse or lTinker et al.l (120081) A = 200 mass functions (see the right panel of Figure |6a]l. 

Changing the overdensity threshold in the spherical overdensity halo definition has an intuitive effect on the linear overdensity 
of regions that collapse to form halos. A higher threshold overdensity, such as a A = 600 halo mass definition (Figurel6b]i. results 
in a higher characteristic linear overdensity for regions that collapse to form halos. Model fits to the mean overdensity with 
smoothing scale show that the mean overdensity increases roughly linearly, and approaches 6 ^ 6c as a ^ 0. The abundance of 
halos is correspondingly suppressed, with the model barrier fits producing very similar excursio n set first-crossing distributions 
that lie below the spherical collapse mass function. The simulated A = 600 mass function from Tinker et"aLl (l2008l) has a lower 
abundance than the lower threshold A = 200 mass function, but the excursion set mass functions are significantly lower than the 
simulated A = 600 mass function (Figure|6b] right panel). In this case, increasing the overdensity threshold does not improve the 
performance of the excursion set ansatz- Similarly, using a lower threshold such as a A = 100 halo mass definition (Figure l6c]i. 
does not bring the excursion set and simulated mass functions into agreement. The lower overdensity threshold decreases the 
characteristic linear overdensity of regions that collapse to form halos and increases the abundance of regions that can collapse. 
However, the mean overdensity does not decrease significantly below the spherical collapse overdensity and the excursion set 
mass functions only begin to roughly match abundance predicted by the spherical collapse mass function. The simulated A = 100 
mass function has increased the abundance of halos relative to the A = 200 mass function, so the disagreement between the 
excursion set mass function and the simulated mass function still remains for this lowered overdensity threshold (Figure|6c] right 
panel). 

Altering the mass definition from spherical overdensity to Friends-of-Friends (FOF) halos does not improve the agreement 
between the simulated and excursion set mass function. The original motivation presented by lSheth et alJ (1200 ll) for modifying 
the ellipsoidal collapse barrier to the lower limiting value of (5smt ^ ^asMT'^e as cr ^ was the use of a FOF definition in 
identifying halos in the GIF simulation. Figure |6d] demonstrates explicitly that the characteristic linear overdensity of regions 
that collapse to form FOF halos does not follow Equation |20l the mean overdensity is similar to that found for A = 200 halos. 
Similarly, the excursion set mass functions calculated from the model fits to the mean overdensity of regions that collapse to form 
FOF halos do not match the simulated FOF halo mass function. 

We could also redefine the mass from the A = 200 halo definition to include all particles within a modified "vi rial" radius 
R = IR yno- This halo mass definition is intended as an analogy to the "static" mass halo definition proposed by Cuesta et"an 
( I2OO8I) . who found that halos defined by regions with zero mean radial velocity (with a size of approximated 2 x /?vir, see their 
Figure 14) displayed an abundance similar to the Pres s-Schechter spherica l collapse mass function at z = 0. This halo definition 
results in halo abundance roughly twice that found by lCuesta et al.l (120081) at fixed mass because the effective Rni/Rioo ratio is 
mass-dependent, but provides a useful example of a mass definition that incorporates very large regions into single halos. Figure 
|6i]shows that the linear overdensity of regions that collapse to form halos defined in this manner is typically low, since the typical 
overdensity at 2/?20o is quite low, and decreases below the spherical collapse barrier at small cr. Hence, we only fit the mean 
overdens ity w ith a linear barrier model that allows for 6 < 6c and do not report the best fit parameters for the model defined by 
Equation lDll As is clear from the right panel of Figure|6i] the excursion set mass function calculated for the linear barrier model 
fit for the R = 2/?20() halos does not recover the simulated halo mass function and, therefore, this halo mass definition does not 
improve the success of the excursion set ansatz- 

Lastly, as in the discussion in ^ one might wonder whether the disagreement between the simulation and excursion set mass 
functions simply involves a normalization issue. In fact, the normalization and shape of each of the excursion set mass functions 
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TABLE 1 

Best Fit Barrier Model Parameters 



Mass Definition 


/3 


7 


A 


b 


A= 100 


0.396 


1.242 


0.442 


1.441 


A = 200 


0.411 


0.809 


0.487 


1.506 


A = 600 


0.543 


0.496 


0.576 


1.660 


FOF 


0.363 


0.890 


0.458 


1.493 
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Fig. 6a. — Smoothed linear overdensity S, extrapolated to z = 0, as a function of smoothing scale a(M) for regions that collapse to form A = 200 halos by 
the present epoch (left panel). The circles coiTespond to the mean overdensities, while the errorbars indicate the halo-to-halo scatter. The en'or on the mean is 
significantly smaller than the scatter in all cases. Shown for comparison is the spherical collapse barrier (5c, solid gray line). Solid red line shows a fit of the 
functional form of the ellipsoidal collapse barrier to the simulation result s, while the da s hed g reen line shows a simple linear fit. The right panel shows excursion 
set mass functions for each model barrier calculated using the method of lZhang & Hull 420061) with sharp k filter and compared with the spherical collapse (gray 
line, right panel) and . Tinker et al.. t.200Sl A = 200 (blue line, right panel) mass functions. 



differ from the mass function constructed from the simulated halo abundance. Figure |7] shows fsimWc)/.fnti'^c)^ the ratio of the 
first-crossing distribution /sim(t'e) determined from the simulations to /fit(j^c), the first-crossing distribution calculated using the 
linear barrier model fits to the overdensity distribution, for each of the halo mass definitions considered in this paper (A = 100, 
A = 200, and A = 600 spherical overdensities, FOF halos, and halos with masses determined by the particle distribution within 
a radius R = 27? 2 00 of the halo center-of-mass). For the A = 200 and A = 600 halo definitions. Figure [T] shows the ratio of the 
best fit analvtical lTinker et al.l (|200 8) and excursion set first-crossing distributions. For the other halo definitions, the binned first- 
crossing distribution determined by the simulated halo mass function is divided by the excursion set first-crossing distribution 
at the appropriate 1/^ value. Uncertainty estimates for fsimit^c) / fmi^'c) for these halo mass definitions (the error bars in Figure |7]i 
have the same fractional error as for the first-crossing distributions plotted in Figures I6cll6el For each mass definition, /sim/./fit 
varies with the peak height and demonstrates that the disagreement between the simulated and excursion set first-crossing 
distributions does not owe simply to their relative normalization. 
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o-(M) V = d /o-(M) 

c c 

Fig. 6b. — Same as for Figure|6a| but for halos defined with a A = 600 spherical overdensity criterion. Increasing the overdensity tlireshold in the halo mass 
definition does not improve the agreement between the excursion set and simulated halo mass function. 




Fig. 6c. — Same as for Figure|6a] but for halos defined with a A = 100 spherical overdensity criterion. Here, the simulated halo mass function is measured 
directly from only the LIOOOW simulation (points and error bars). 
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o-(M) V = d /o-(M) 

c c 

Fig. 6d. — Same as for Figure|6c] but for halos defined with a Friends-of-Friends (FOF) algorithm using a hnking length of b = 0.2. The chara cteristic linear 
overdensity of regions that collapse to form FOF halos does not follow the modified ellipsoidal collapse barrier presented by dSheth et"aI]|2001L purple dotted 
line). Changing the halo definition from spherical overdensity to FOF halos does not improve the agreement between the excursion set and simulated halo mass 
function. The simulated halo mass function is measured directly from only the LIOOOW simulation (points and en'or bars). 




ff(M) V = d /<7(M) 

c c 

Fig. 6e. — Same as for Figure l6cl but for A = 200 halos with masses rescaled to include all particles within a radius R = 2S200- In this halo definition, small 
A = 200 halos with centers-of-mass that reside within a distance R of larger halos are incorporated into the larger system. The mean of the linear overdensity 
distribution in this case lies below the spherical collapse bairier at small cr, so only the excursion set mass function for the linear bairier (which allows S < 5^) 
is compared with the simulated mass function (right panel). Changing the halo mass definition to increase the region incorporated into halos identified by a 
A = 200 spherical overdensity criterion does not improve the agreement between the excursion set and simulated halo mass function (points and eiTor bars). 
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Fig. 7. — Ratio /sim(i'i:)//fi(i^c) of the simulated first-c rossing distribution /sip, to tlie excursion set first-crossing distribution /fjt calculated from a linear 
fit to the overdensity distribution plotted in Figures l6al6el for a variety of halo mass definitions, plotted as a function of peak height Uc = Sc/u(M). For the 
A = 200 (orange fine) and A = 600 (red line) spherical overdensity halo definitions, the analytical fits fromlTinker et alj (12008) ai'e compared with the calculated 
first-crossing distributions. For A = 100 spherical overdensity halos (purple points), Fiiends-of-Friends halos (blue points), and halos with masses defined by the 
particle content within a radius R = IRiqo (green points), the binned first-crossing distributions from the LIOOOW simulation halo catalogue were compared with 
the excursion set results at the appropriate Uc for each mass bin. For the se halo mass definitions, the uncertainty estimates in fsim/ffit reflect the same fractional 
uncertainty in the simulated first-crossing distributions shown in Figures [6c|6e| For every halo mass definition, that /sim/Zflt is not constant with Uc demonstrates 
that the disagreement between the simulated and excursion set first-crossing distributions differ in shape and not simply in relative normalization. 



